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We characterize gapless edge modes in translation invariant topological insulators. We show that 
the edge mode spectrum is a continuous deformation of the spectrum of a certain gluing function 
defining the occupied state bundle over the Brillouin zone (BZ). Topologically non-trivial gluing 
functions, corresponding to non-trivial bundles, then yield edge modes exhibiting spectral flow. We 
illustrate our results for the case of chiral edge states in two dimensional Chern insulators, as well 
as helical edges in quantum spin Hall states. 
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The study of topological phases of matter has been an 
exciting field of research since the discovery of the inte- 
ger quantum Hall effect (IQHE) in the 1980's. Recently, 
the discovery of the quantum spin Hall effect (QSHE) 
[H-Q and three dimensional topological insulators 0d3] 
has shown that interesting physics occurs even in these 
simple band models. The common feature uniting these 
materials is a topological "twisting" of the band struc- 
ture over the BZ; stated mathematically, the invariant is 
the K-theory class of the occupied state vector bundle 
over the BZ [9|, |l0(. Another striking characteristic of 
these materials is gapless edge modes. In specific cases, 
like the IQHE or even QSHE their existence is guaran- 
teed by various arguments 11-lH but it is natural to 
ask whether there is a more direct connection between 
bulk invariants and gapless edge modes. Certainly, there 
are cases of non-trivial topology, such as stacked IQHE 
planes, where edges cut parallel to the planes can be 
gapped. Given that non-trivial topology by itself does 
not imply gapless edge modes, when do protected gap- 
less edge modes occur? 

In this Letter, we give a necessary and sufficient con- 
dition for protected gapless edge modes in terms of topo- 
logical data that define the occupied state bundle. We 
study a planar edge perpendicular to a crystal direc- 
tion, giving a preferred splitting of the BZ torus T d into 
d — 1 "parallel" directions and a "perpendicular" direc- 
tion: T d = T d ~ x x S 1 . With respect to this splitting, 
one may regard the occupied state bundle as a bundle on 
T d_1 , extended trivially to T^ 1 x [0, 2n], and then glued 
along the boundary (d — l)-tori at and 2ir. This gluing 
is encoded in a function U g : T d ~ 1 — > U(N), where N is 
the number of occupied bands. The non-Abelian Berry 
connection gives a natural way to locally straighten the 
fibers in the perpendicular direction, and U g can then 
be defined via parallel transport along the perpendicular 
S 1 : 



where A± is the perpendicular component of the Berry 
connection, and the exponential is path-ordcrcd. The 
gluing function U g has a simple interpretation in terms of 
localized Wannier functions [16| : being the exponential of 
the perpendicular Berry-covariant derivative in momen- 
tum space, U g is, in real space, the exponential of the per- 
pendicular position operator, projected into the occupied 
bands. The corresponding eigenfunctions are localized 
Wannier states |Xn( r ~ ' a ))i where a is the unit cell spac- 
ing. Suppressing parallel components, the \xn{ r — lei)) 
are related to Bloch functions \u n (k)) via a generalized 
Fourier transform 



\Xn(r-la)) = J 



dk ± e lk ^ r ~ l ^\u n (k)). 



(2) 



2 71 



U g (k\\) = exp I i / A±(k\\,kx)dk± 



(1) 



The eigenvalues of U g are then of the form exp(27rz</>„), 
where 1 < n < N and <f> n , defined modulo 1, gives real 
space position of the center of the Wannier function Xn 
modulo the unit cell. Our central claim is that the spec- 
trum {<f> n (k\\) + I}, where / ranges over the integers, can 
be continuously deformed into the edge mode spectrum. 
Thus the two spectra have the same topology, so spectral 
flow of edge modes is equivalent to spectral flow of the 
Wannier centers. Because the latter is a property only of 
the bulk, this connection yields a criterion for the exis- 
tence of protected gapless edge modes purely in terms of 
the bulk band structure. After deriving the general re- 
sult, we illustrate it by example for Chern insulators and 
QSH systems; in general, it can be applied to systems in 
any dimension and symmetry class. 

When constructing the edge, we have a choice of 
boundary conditions, but all such choices are deformable 
to each other and hence, in the presence of a bulk gap, 
have the same topology. Furthermore, in the bulk the 
gapped Hamiltonian can be continuously deformed to one 
which is spectrally flat, without closing the bulk gap and 
hence without changing the topology of the edge states. 
Thus, we start with a spectrally flat tight binding Hamil- 
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tonian H = 1 — 2P defined on an infinite lattice Z , where 



P 



a/3 



B.Z. 



d k J(r-s) p af) 

{2ttY k ' 



(3) 



r, s £ Z d are discrete spatial lattice indices, k £ T d the 
dual momentum, and a, j3 = 1, . . . , N the band indices. 
Pfc satisfies P| = Pk , and is thus a projector denning the 
occupied states at momentum k. We choose a convenient 
way to introduce an edge to the system by defining 



PV (x)P + (1 - P). 



Here 



i 1, for xx < 
W = { _j for ^ > 



(4) 
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and xj_ is the component of x perpendicular to the edge. 
H e has exponentially decaying hopping because P = (1 — 
i?)/2 does, and Vo is local. To the right, it approaches 
the spectrally flat form H lig h t = 1 — 2P exponentially 
fast, while to the left it approaches the trivial atomic 
limit with negative chemical potential, H\ e f t = 1. Thus 
it models a boundary between topological insulator and 
vacuum. Since P and 1 — P act in orthogonal blocks, the 
entire spectrum of H e , aside from the trivial eigenvalue 
1, is contained in PVq(x)P. 

To analyze the spectrum of PVo(x)P, we first note that 
fcii , the momentum parallel to the cut, is conserved, so 
that we can study each fc|| sector separately. We first 
prove that the spectrum of PVo(x)P, in each kn sector, 
is discrete. These discrete levels correspond to edge sub- 
bands; in contrast to a generic gapped Hamiltonian, the 
special one defined in (j4]) has no band continuum, but 
rather only accumulation points at ±1, corresponding to 
edge sub-bands penetrating further and further into the 
bulk. Demonstrating this for fixed ku is effectively a one 
dimensional problem, so to avoid cumbersome notation 
we simply drop feii and assume we have a one dimensional 
system in the following proof. 

It is useful to view P as a matrix with respect to the 
partition of the Hilbert space to the left and right of the 
cut: 



P 



Pll Prl 
Plr Prr 



(6) 



Plr, Prl are compact operators, since their matrix ele- 
ments decay exponentially, implying that Pll + Prr = 
P — Plr, — Prl has a discrete spectrum with possible 
accumulation points at 0,1. Since Pll and Prr act in 
orthogonal subspaces, they have discrete spectra as well, 
as do Hll/rr = 1 — ^Pll/rr- The latter have accumu- 
lation points at ±1. Elementary computation now gives 



PV (x)P-(l-P)V (x)(l-P) = 



-H LL 
H RR 



(7) 



Since PV (x)P and (1- P)V (x)(l~ P) act in orthogonal 
subspaces, both must have discrete spectra with possible 
accumulation points at ±1, as desired. We note here 
that the spectrum of ffsnis directly related to the en- 
tanglement spectrum 17H25| via Peschel's construction 
24l | . Under mild non-degeneracy assumptions, namely 
the lack of bulk occupied states localized exactly on one 
side of the cut, the spectra of PVq(x)P and Hrr are 
actually the same. 

We now deform the spectrum of PVo(x)P to that of 
(l/27ri) log U g . This deformation will be continuous in 
fc|| and preserve discreteness at each fixed fc||, thus pre- 
serving topological properties as well. For convenience 
we again fix k\\ to work with a one dimensional system 
in the following argument. Recall that (l/2-7r) log C/g has 
spectrum {(/>„ + I}, where < 4> n < 1 and I £ Z. Wc 
first re-interpret this as the spectrum of the Berry gauge- 
covariant derivative — iV b = —idk ± + ^4j_ acting on sec- 
tions of the occupied state bundle. Indeed, defining a 
covariantly-constant frame \u n (k±)} we see, using ([1}, 
that |u„(27r)) = ([/g)™|u rn (0)); we can furthermore en- 
sure that Ug is diagonal in the latter expression. Then 
the eigenfunctions of — iV b are localized Wannier states 
Xn( r — la) = J dk± expifcj_(</>„ + I) \u n (kj_)) which have 
eigenvalues 0„ + I, as desired. Note also that, as an op- 
erator on Hilbert space, — i\7 B = P(—idk ± )P = Px±P. 
Thus, all we have to do is deform the operator PVq(x)P 
to Pxj_P. 

Define 



V t {x) = 



-x±, 

-sgn(xj_)/(l 



for < 1/(1 -t), 
t) for \x±\ > 1/(1 - t). 



(8) 

Then Vt(x) interpolates between Vq(x) and x±, and 
PVt{x)P is the desired deformation. Indeed, for t < 1, 
PVt{x)P is a finite rank perturbation of (1 — t)^ 1 PVqP, 
acting only at \xj_\ < 1/(1 — t), and thus cannot produce 
a continuous spectrum. For t = 1 a slight subtlety arises: 
Px±P is not a bounded operator. To make a rigorous 
statement we define a bounded continuous function h(y) 
such that h(y) — y on [— W, W] and h(y) = sgn(y)W for 
|x| > W, for some large W. Then one can check that the 
spectrum of the uniformly bounded family of operators 
h(PVtP) is jointly continuous in t £ [0,1] and k\\. This 
is just the statement that as one goes from PVq(x)P to 
PV\(x)P, the spectrum evolves uniformly continuously 
in any finite window [— W, W]. 

We have thus proven that the spectrum of the loga- 
rithm of Ug and the edge spectrum have the same topol- 
ogy. We now illustrate our result for Chern insulators 
and time reversal invariant systems in two dimensions. 

Chern insulators: Time reversal breaking insulators in 
two dimensions are characterized by an integer Chern 
number (TKNN invariant). The simplest realization of 
nonzero Chern number has one filled band, corresponding 
to a line bundle over the T 2 BZ. The transition function 
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discussed above is a map from S 1 C T 2 to U(l), and the 
element of tti(S 1 , U(l)) = Z which it defines is equal to 
v. The spectrum of its logarithm has v protected modes 
spiraling up, just as the edge spectrum, consistent with 
our picture - see figure [JJ 

More quantitatively, we can see spectral flow directly 
by evaluating [2(| 



u = 2niT r (P[[P ) U x ],[P,U y }} 



(9) 



where now H x ,y are projections onto the positive x and 
y half-planes respectively. The topological nature of v 
gives us a lot of freedom to deform the expression ([9]). 
In particular, with translation symmetry one can deform 
the Hx,y m to linear functions X, Y (replacing sums with 
integrals), obtaining the usual expression for v as the 
integral of the Berry curvature over the BZ j26[. Now 
take an entanglement cut along y (so that x = xj_ and 
y = .T||). Then it is most useful to deform only Tl y : 



v = / dk y TrC(k y ) 



(10) 



where C{k y ) = P ky [[P ky ,U x ],dP ky /dk y }. We evaluate 
(ITUf in the basis \ip n (k y )) of eigenstates of P ky Vo(x)P ky : 



Pk y V (x)P ky \ip n ) = \ n \i> n ). 



(11) 



Using (J7]) and the fact that -^rH^Vi) = gf-IVVi) is or- 
thogonal to |"0ri) we obtain 







WC(*»)IW = (V'nin.lVn) = — 



d fl + K 



dk q , 



Thus IjlOp is the integral of a total derivative, resulting 
in the expression v = ^2 n (^n(2TT) — A n (0))/2. Now, 
A„(27r) = A n+ / £ (0) for some k which is easily seen to 
be independent of n (the case of degenerate A„ must be 
handled carefully, but all level crossings can be appropri- 
ately resolved), k is simply the number of edge modes 
that have to cross any given Fermi level. The sum then 
telescopes and we obtain v = fc, as desired. 

Continuum Chern insulators: [27, 28 1 As a check of the 
above results, one would like to have a direct calculation 
of the spectrum of Hrr = 1 — 2Prr- This is tedious 
for a lattice Chern insulator, but easily doable in the 
continuum, i.e. for the IQHE. One problem is that the 
continuum BZ is not a finite torus; nevertheless, we will 
satisfy ourselves with demonstrating spectral flow as fen 
ranges over the momentum scale corresponding to the 
inverse magnetic length i~ x . We will work on a semi- 
infinite cylinder of large radius; this only introduces a 
fine discretization of the domain of the plots. 

The bulk Hamiltonian H = 1 — 2P gives energy —1 to 
the v lowest occupied Landau levels, and energy 1 to all 
higher levels. Its eigenvalues are labeled by fcy and a Lan- 
dau level number, and the eigenstates are localized har- 
monic oscillator functions in the xj_ direction. The trun- 
cated Hamiltonian Hrr, written in this basis, contains 




FIG. 1: Spectrum of gluing function (top) and entanglement 
- or edge - spectrum (bottom) for a) a Chern insulator with 
one occupied band and v — 1 and b) a minimal quantum spin 
Hall system with nontrivial Z2 invariant. 



Landau level mixing terms. Letting n, m = 0, . . . , v — 1 
denote occupied Landau levels, (Hjm)^ n is equal to 



i\m\ 



due" 2 (c£e-« 2 ) (cCe^ 2 ). 



(12) 

Here uq = k\\£ and u = x±/£ + uq- For example 

For example, for v = 2 the eigenvalues of (Hrr) 7 ^ are 
plotted in Fig. [5J we see that both branches cross any 
given value of A, displaying spectral flow. 




FIG. 2: Entanglement spectrum flow Eq. (??) for v — 2 
Landau levels. 

Quantum spin Hall systems: QSH insulators are char- 
acterized by a nonzero Z2 invariant. With a choice of 
gauge one can express it in terms of Pfaffians [l| , or as 
an obstruction to trivializing the occupied state bundle 
in a way compatible with T jl6| . These are equivalent to 
the gluing function constructed above being non-trivial. 
Indeed, since the spectrum of the gluing function is given 
by the positions of the Wannier centers, a non-trivial glu- 
ing function is equivalent to Kramers pair-switching spec- 
tral flow of the time reversal polarization, an equivalent 
formulation of the Z2 invariant [lij . 

Because T symmetry relates H (fc) to H(—k), the spec- 
trum of U g (ku) is the same as that of U g {— fcii), and be- 
cause T 2 = — 1 requires an even number of bands, the 
gluing function can be thought of as a map U g : [0, n] — > 
U(2N). The extra constraint TVgT- 1 = Ug 1 at the T 
symmetric endpoints fcii = 0, tt forces the spectra at those 
points to form Kramers pairs. The non-trivial class now 
switches elements in these Kramers pairs as one moves 
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from to 7r (formally the classification is given by the 
relative homotopy group m(U(2N), U(N)) = Z 2 ). 

The key point now is that the deformation between 
the spectrum of ^ log U g and the edge spectrum re- 
spects T throughout, so Kramers partner-switching in 
the edge spectrum is an equally good criterion for QSH 
topological order. The simplest example consists of two 
time-reversed filled bands, with the two eigenvalues of 
the gluing function splitting up as one moves away from 
k y = 0, traversing by ir in opposite directions and re- 
joining at k y = 7r; the helical edge spectrum shares this 
property, consistent with our result (Fig. Q}. 

Generalizations: 

We related only one particular gluing function to edge 
modes. By considering the other d — 1 possible cuts, we 
gain more information. In particular, if the edge modes 
for all d cuts are gapped, our construction shows that 
the insulator is topologically trivial, i.e. can be adiabati- 
cally connected to a trivial insulator in the atomic limit. 
Furthermore, we can perform similar arguments in other 
dimensions and symmetry classes, e.g. 3d topological in- 
sulators pij . 

Another generalization concerns our formula for the 
Chern invariant in terms of the entanglement spectrum. 
We can in fact apply it to any mixed state p that is trans- 
lationally invariant in y. What does such an invariant 
represent? 

To answer this we construct the purification of our den- 
sity matrix. Consider p = Z~ x e~ ^ H rr / *<"-*0»-' _ We now 
double the number of sites and consider the Hamiltonian 



HAW — J2rar'a( 



I) Q ^ t q, R . + (1 

2 / rr' ar * pr' ' V 2 



+ v/n(l - nfj, ^L^r* + V /^T 



where n = j^pr is the Fcrmi-Dirac operator and are 
the fermi operators on the auxiliary system. Writing in 
block form the single particle Hamiltonian of Haw, it is: 



v /n(l - n) 



y/n(l - n) 



It is easy to check that hwA is of the form P — |, P a 
projector. Hence, Haw is spectrally fiat, describing two 
bands of energies —1/2 and 1/2. The ground state is 
the purification of our state p, also known as the "Araki- 
Wyss" representation [3(J. For a system with decaying 
two point function, one can check that n has decaying 
matrix elements. Thus the Haw above describes a bi- 
layer system coupled in the bulk, rather than through an 
edge. Indeed for a generic thermal state the entangle- 
ment entropy associated with H is extensive. However if 
H is coming from tracing half of a system at zero temper- 
ature, the coupling terms, proportional to \/n(l — n) will 
vanish exponentially inside the bulk, since n(l — n) = 
whenever n = or n = 1, and we are left with a bilayer 
which is only coupled close to the boundary. Thus, the 



entanglement Chern number represents the Chern invari- 
ant of the associated bilayer Araki-Wyss system. 

Conclusions: We have shown that the edge modes of 
a topological insulator are a continuous deformation of 
the spectrum of a gluing function defining the occupied 
state bundle over the BZ, giving a necessary and sufficient 
condition for protected gapless edge modes. 

We thank B. A. Bcrncvig, M. Frccdman, and espe- 
cially A. Turner for useful discussions. IK acknowledges 
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